General Case: Eigenvalues and Eigenvectors 



We are now ready to tackle the general case of a linear 2x2 system: 

x — ax + by 
y = cx + dy , 

where the a, b, c, d are constants. We will be following exactly the strategy 
that we laid out in the previous note. These are key concepts for the rest of 
the unit, and you should take the time to absorb them. 

We want to learn to write the system efficiently in matrix form. So, 
throughout the derivation, we will give the expanded matrix form of our 
manipulations on the left, and the abridged form on the right. For example, 
our system is: 

*)-(:!)(;) - 

Following the method in the previous note, we look for solutions to our 
system having the form 

where a\, ai, and A are unknown constants. We substitute this into the 
system (1) to determine these unknown constants. Since D(ae M ) = Aae M , 
we arrive at 

( Ul ) =e At ( U b A( Ul ) * A**a = e kt A ( Ul 
We can cancel the factor e M from both sides, getting 

A / a x \ = / a b \ / a x \ ^ \a = A ( Ul 

As it stands, we cannot combine the two sides by subtraction, since the 
scalar A cannot be subtracted from the square matrix on the right. As in 
the previously worked example, the trick is to replace the scalar A by the 
diagonal matrix XL This gives 

i / 1 ■ 11 '- - A /a = A 1 
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If we now proceed as we did in the example, subtracting the left side from 
the right one and using the distributive law for matrix addition and multi- 
plication, we get a 2 x 2 homogeneous linear system of equations: 



Written out without using matrices, the equations are 

(a — A)fli + ba-i = 
ca\ + (d — A)fl2 = 0. 



/ 



(2) 



According to the theorem on square homogeneous systems this system has 
a non-zero solution for the a's if and only if the determinant of the coeffi- 
cients is zero, i.e., 



a — A b 
c d — A 



= 



\A-\I\ = 0. 



Evaluating the determinant we get a quadratic equation in A: 

A 2 - (a + d)A + {ad -be) = . 



Definition. This is called the characteristic equation of the matrix 




and if often denoted Pa (A). Its roots Ai and A2 are called the eigenvalues 
or characteristic values of the matrix A. 

Remark. In calculating the characteristic equation notice that 
ad — be = det A a + d = tr A. 

Using this, the characteristic equation for a 2 x 2 matrix A can be written 

as 

A 2 — (tr A) A + det A = 0. 

In this form, the characteristic equation of A can be written down by inspec- 
tion; you don't need the intermediate step of writing down I A — AI| = 0. 

Remark. Abridged vs. expanded notation 

In the manipulations above, the matrix notation on the right is compact 



2 



General Case: Eigenvalues and Eigenvectors 



OCW 18.03SC 



to write, which makes the derivation look simpler. On the other hand, 
its chief disadvantage for beginners is that it is very compressed. Practice 
writing the sequence of matrix equations so you get some skill in using the 
notation. Until you acquire some confidence, keep referring to the written- 
out form on the left, so you are sure you understand what the abridged 
form is actually saying. 

There are now various cases to consider, according to whether the eigen- 
values of the matrix A are: 

1. two distinct real numbers, 

2. a single repeated real number, 

3. a pair of conjugate complex numbers. 

We begin with the first case: for the rest of this note, the eigenvalues are 
two distinct real numbers K\ and A2. 

1. Real distinct eigenvalues 

To complete our work, we have to find the solutions to the system (2) 
corresponding to the eigenvalues k\ and A2. Formally, the systems become 



The solutions to these two systems are column vectors, for which we will 
typically use v. 

Definition. The respective solutions a = vi and a = V2 to the systems 
(3) are called eigenvectors (or characteristic vectors) corresponding to the 
eigenvalues Ai and A2. 

Remarks. 1. If the work has been done correctly, in each of the two sys- 
tems in (3), the two equations will be dependent, i.e., one will be a constant 
multiple of the other. Why? The two values of A have been selected so that 
in each case the coefficient determinant A — AI will be zero, which means 
the equations will be dependent. 

2. The solution v is determined only up to an arbitrary non-zero constant 
factor: if v is an eigenvector for A, then so it cv, for any real constant c; 
because of this, the line through v is sometimes called an eigenline. A 
convenient way of finding the eigenvector v is to assign the value 1 to one 
of the di, then use the equation to solve for the corresponding value of the 
other a.{. (First try a\ = 1; if that does not work, then ai = 1 will.) 

Once the eigenvalues and their corresponding eigenvectors have been 



(a — Ai)fli + bai — 
ca\ + (d — Ai)fl2 = 



(a — X-i)a\ + bai = 

Cfli + (d — \2) a 2 — 



(3) 
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found, we have two independent solutions to the system (1). They are 



Definition. In science and engineering applications, these are usually 
called the normal modes. 

Using superposition, the general solution to the system (1) is 



Remarks. 1. The normal nodes often have physical interpretations; this 
means that it is sometimes possible to find them just by inspection of the 
physical problem. 

2. In the compact notation, the definitions and derivations are valid for 
square systems of any size. Thus, for instance, you know how to solve a 
3 times?> system, if its eigenvalues turn out to be real and distinct. We won't 
consider any such systems in these notes, though. 

We will apply these techniques to a worked example in the next note in 
this session titled Worked Example: Distinct Real Roots. 



xi(0 = e Al S, 




X = C\X\ + C2X2 = Ci6 Alf Vi + C2le k2t \2- 
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